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ABSTRACT 


A crucial part of aircraft design is ensuring that the required margin for flutter is satisfied. A 
trustworthy flutter analysis, which begins by possessing an accurate dynamics model, is necessary 
for this task. Traditionally, a model was updated manually by fine tuning specific stiffness parameters 
until the analytical results matched test data. This is a time consuming iterative process. The NASA 
Dryden Flight Research Center has developed a mode matching code to execute this process in a 
more efficient manner. Recently, this code was implemented in the F- 1 5B/Quiet Spike™ (Gulfstream 
Aerospace Corporation, Savannah, Georgia) model update. A build-up approach requiring several 
ground vibration test configurations and a series of model updates was implemented to determine 
the connection stiffness between aircraft and test article. The mode matching code successfully 
updated various models for the F-15B/Quiet Spike™ project to within 1 percent error in frequency 
and the modal assurance criteria values ranged from 88.51-99.42 percent. 

NOMENCLATURE 


CAD 

computer aided design 

CG 

center of gravity, in 

CONM2 

concentrated mass element 

F 

temporary variable 

FE 

Finite Element 

G 

target value for optimization 

GAC 

Gulfstream Aerospace Corporation 

GVT 

ground vibration test 

i 

index term 


x-axis, lb-in 2 
x-axis, lb-in 2 , target 
y-axis, lb-in 2 
y-axis, lb-in 2 , target 
z-axis, lb-in 2 
z-axis, lb-in 2 , target 
y-axis, lb-in 2 
y-axis, lb-in 2 , target 
z-axis, lb-in 2 
z-axis, lb-in 2 , target 
x-axis, lb-in 2 
x-axis, lb-in 2 , target 

j second index term 


^XX 

moment of inertia 

^XXG 

moment of inertia 

I yy 

moment of inertia 

^YYG 

moment of inertia 

^ZZ 

moment of inertia 

^ZZG 

moment of inertia 

^XY 

moment of inertia 

^XYG 

moment of inertia 

Iyz 

moment of inertia 

^YZG 

moment of inertia 

^ZX 

moment of inertia 

^ZXG 

moment of inertia 


about the x-axis with respect to the 
about the x-axis with respect to the 
about the y-axis with respect to the 
about the y-axis with respect to the 
about the z-axis with respect to the 
about the z-axis with respect to the 
about the x-axis with respect to the 
about the x-axis with respect to the 
about the y-axis with respect to the 
about the y-axis with respect to the 
about the z-axis with respect to the 
about the z-axis with respect to the 


J objective function 

K stiffness matrix 

M mass matrix 

MAC Modal Assurance Criteria 

m number of degrees of freedom 

n number of modes that will be correlated in eq. (2) 

NaN not a number 

QS Quiet Spike™ 

R temporary variable 

T transformation matrix 

x-CG center of gravity, x-axis 

y-CG center of gravity, y-axis 

|3 scaling factor 

O analytical Eigen matrix 

to. test frequencies 

Q. analytical frequencies 

V temporary variable 


INTRODUCTION 

When a new aircraft is designed or a modification is made to an existing aircraft, the aeroelastic 
properties of the aircraft need to be examined to ensure it is flightworthy. Evaluating the aeroelastic 
properties of a new or modified aircraft can include performing a variety of analyses, such as modal 
and flutter analyses. To produce accurate results from these analyses, it is imperative to work with 
finite element (FE) models that have been validated by or correlated to ground vibration test (GVT) 
data. Updating an analytical model using measured data is a challenge in the area of structural 
dynamics. The analytical model update process encompasses a series of optimizations that match 
analytical frequencies and mode shapes to the measured modal characteristics of a structure. In the 
past, the method used to update a model to match testing data has been ‘trial and error.’ This is an 
inefficient method; running a modal analysis, comparing the analytical results to the GVT data, 
manually modifying one or more structural parameters [mass, center of gravity (CG), inertia, area, 
etc.], rerunning the analysis, and comparing the new analytical modal characteristics to the GVT 
modal data make up a time-consuming process. If the match is close enough (close enough defined 
by the analyst’s updating requirements), then the updating process is finished. If the match does not 
meet updating requirements, then the parameters are changed again and the process is repeated. 
Clearly, this manual optimization process is highly inefficient for large FE models and/or a large 
number of structural parameters. 

The NASA Dryden Flight Research Center (DFRC) has developed, in-house, a mode 
matching code that automates the above-mentioned optimization process. Moreover, it facilitates 


the development of reduced-order models that simplify analyses and therefore, reduce computation 
time. The DFRC in-house mode matching code reads mode shapes and frequencies that were 
measured from a GVT to create the target model. It also reads the current analytical model as the 
‘initial guess,’ as well as the variable parameters and their upper and lower limits as assigned by 
the analyst. It performs a modal analysis on this initial model and modifies its variable properties 
to create an updated model that has similar mode shapes and frequencies to those of the target 
model. The mode matching code outputs frequencies and modal assurance criteria (MAC) values 
that allow for the quantified comparison of the updated model versus the target model. 

The mode matching code assists in determining the airworthiness of the aircraft flown 
at DFRC. Recently, this code was implemented on a project flown on the DFRC F-15B (The 
Boeing Company, Chicago, Illinois) supersonic flight-testing platform, the Gulfstream Aerospace 
Corporation (Savannah, Georgia) (GAC) Quiet Spike™ (QS). This code allowed for a build-up 
approach that was necessary because the flight test article itself was unavailable. Researchers were 
able to make progress in the analyses even without the flight hardware. 

The objective of this work was to use the mode matching code to determine the connection 
stiffness on the airplane radar bulkhead where the QS and the F-15B airplane mate, such that when 
the FE model of the flight QS was validated, it could be readily mated to the airplane FE model. 
This approach was established to enable the project to progress the model development and flutter 
predictions without the actual QS flight test hardware. 

MATHEMATICAL BACKGROUND 

Three optimization phases are used to refine a model with the DFRC mode matching 
code. The mass properties are modified to match measured mass properties, the mass matrix is 
orthogonalized, and then the natural frequencies and mode shapes are matched. Design variables 
for the optimization can include structural sizing information (thickness, cross-sectional area, area 
moment of inertia, torsional constant, etc.), point properties (lumped mass, spring constant, etc.), 
and material properties (Young’s modulus, etc.). 

Phase 1: Mass properties 

Matching mass properties by using lumped masses (CONM2)(ref. 1 ) reduces the computational 
effort required to match stiffness properties and ensures that a physically meaningful and accurate 
solution is achieved. In this first optimization phase, 10 objective functions are minimized to 
match analytical to measured mass properties. The mass properties that are updated include 
total mass, center of gravity location, and mass moments of inertia. As an objective function is 
minimized, a constraint function is implemented to prevent the updated variable from changing 
in the minimization of the following objective function. The constraint functions limit the amount 
of change allowed in the previously optimized variables. The objective functions and constraints 
applied in this phase are listed in order in table 1 . This optimization phase places the model within 
the vicinity of a physically viable solution (feasible region), facilitating the second and third 
optimization phases. 
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Table 1. Optimization problem statement for phase 1 of model update. 


Statement 

Objective Function 

Constraint Function 

1 

y 1 = ( w-w g ) 2 /w g 2 

Unconstrained 

2 

j, = (x-x a y/x a 2 

l/jl < e 



3 

J , = (Y-Y 0 )W 0 2 

I/I < 8 

l 

i = 1, 2 


4 

j, = (Z-Z a ) ! /Z 0 ’ 

I/I < 8 

1 

i = 1, ••• 

3 

5 

j = i i -i y/i 2 

3 v XX XXG' XXG 

I/I < 8 

1 

i = 1, ••• 

4 

6 

J =(l -I WI 2 

6 v YY YYG' YYG 

I/I < 8 

1 

i = 1, 

5 

7 

J = (I -I ) 2 /l 2 

/ v ZZ ZZG' ZZG 

I/I < 8 

1 

i=l,... 

6 

8 

J =(l -I ) 2 /l 2 

u 8 v XY XYG^ XYG 

I/I < 8 

1 

i= 1 

7 

9 

J =(l -I ) 2 /l 2 

9 v YZ YZG' YZG 

I/I < 8 

1 

i = 1, ••• 

8 

10 

J = (l -I ) 2 /l 2 

10 v ZX ZXG' ZXG 

I/I < 8 

1 

i= 1 

9 


Phase 2: Mass Matrix Orthogonalization 

The orthonormalized mass matrix is shown in eq. (1). 

M = 3>qT T MTO g (!) 

The off-diagonal terms of this orthonormalized matrix are reduced to improve mass 
orthogonality. In eq. (1), the mass matrix, denoted as M, represents the analytical mass matrix, 
while represents the target Eigen matrix that is obtained from the GVT. The target Eigen 
matrix, <D G , remains unchanged during the optimization. The transformation matrix, T, is based on 
Guyan reduction (ref. 2) or improved reduced system (ref. 3). The optimization problem statement 
then becomes: 

Minimize 


£ (Mi /) 2 (2 > 

(Off-diagonal terms of M ) 

Such that I/I < e where i = 1, 2. . . 10 (/. values are given in table 1), and n denotes the number of 
modes that will be correlated in eq. (2). The positive-definiteness characteristic of the lumped mass 
elements also acts as a constraint during this optimization procedure. 

Phase 3: Frequencies and Mode Shapes 
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In this phase, certain modes can be given a higher level of importance by applying a scaling 
factor, (3, in eq. (3) and eq. (5). This scaling factor helps to more closely tune the frequencies 


of these analytical modes of interest after the mode shape vectors have been matched to within 
required accuracy. 

Two optimization methods are used in this phase. In the first method, the objective function 
considered combines an index which identifies normalized errors between test (to.) and analytical 
(Q.) frequencies with a second index that defines the total error associated with the off-diagonal terms 
of the orthonormalized stiffness matrix. The optimization problem statement then becomes: 

Minimize 


J 12 ~ X 




i=l V £2i J 


+ P 2 (Ky) 


(3) 


Such that I/I < e for i = 1, 2... 11. The matrix K is a result of the matrix multiplication shown in 
eq. (4), where the stiffness matrix, K, is the analytical stiffness matrix. 

K = <D ( It T KTO g (4) 

In the second optimization method, the objective function combines the same frequency error 
index used above, with a second index that defines the total error between test and analytical mode 
shapes at specified sensor locations. The optimization problem statement then becomes: 

Minimize 


J l2 ~ X 


n f n. — m 

CO i 


i=l V J i=l 


+ ^S(Oi-<I>iG) 


(5) 


Such that I/I < 8 for i = 1, 2... 11, and m denotes the number of degrees of freedom measured. 
The frequencies to. and mode shapes O. are calculated based on the direct-iterative Eigen solution 
technique for large eigenvalue problems (ref. 4): 


1) Compute K x and M x : 

2) Compute F: 

3) Compute R: 

4) Compute K R and M R : 

5) Compute to. and V: 

6) Compute O: 


K x = T T KT and M x = T T MT 


F = M x 0> g 


KR = I 

K p = R t K t R and M p = R t M t R 

K 1 K 1 

K P V = to 2 M P W 

K 1 K 

0 = R»P 
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For the direct-iterative Eigen-solution technique, the analytical stiffness and mass matrices, K 
and M, are correlated, while the measured mode shapes, 4> G , remain unchanged during the model 
updating process. Notice that the transformation matrix, T, in eq. (1) and eq. (4) is calculated 
using the analytical mass and stiffness matrices M and K. Therefore, the constraint equation 
U„l < 8 is needed for the third optimization phase. 

Current model updating approaches consist of updating only stiffness properties and assume 
that the modeled mass distribution is accurate. This can create erroneous stiffness properties in a 
model. The DFRC mode matching code is capable of updating both mass and stiffness in a model 
by updating mass properties such as total mass, CG, and mass moments of inertia, improving the 
orthogonalization of the mass matrix, and matching frequencies and mode shapes to test data. 
Updating both mass and stiffness properties of a model will provide a higher degree of confidence 
in a model’s resemblance of the actual hardware. 

APPLICATION 

An example of an application of the mode matching code is the F-15B research testbed. 
The DFRC possesses a modified F-15B airplane that is used as a testbed for supersonic flight 
experiments. Traditionally, the FE model of the flight test article is generated and updated by GVT 
data and then mated to the updated F-15B FE model. A GVT is then conducted with the test article 
mated to the airplane. The mated FE model is correlated to this mated GVT data with the going-in 
assumption that, at this point, the only unknown variable is the interface stiffness between the test 
article and the airplane. 

Recently, GAC expressed a desire to test a modified nose, the Quiet Spike™ , on the DFRC 
F-15B (ref. 5) airplane. Because of flight test article availability and a compressed project schedule, 
a non-traditional-testing and model update method had to be implemented. A build-up approach 
was devised in which a Mock QS was designed and fabricated to replace the flight QS in the tests 
and analyses. A detailed FE model of this Mock QS was created. A reduced-order model of this 
detailed FE model was created in the form of a beam model to expedite the model update process. 
The beam model was then matched to data provided by the first GVT using a strongback. This 
beam model was mated to the F-15B model and the interface stiffness was updated by a second 
GVT using the airplane. The flight QS was tested and its model was updated separately at GAC. It 
was then mated to the aircraft model and its modal analysis results were compared to frequencies 
and mode shapes from a final GVT of the flight QS mated to the F-15B airplane. Ultimately, the 
goal of this build-up approach was to completely eliminate the need for a model update of this final 
flight test configuration. 

The flight QS measures 30 ft in length when fully extended and 20 ft when fully retracted. 
Figure 1 shows the GAC flight QS. 
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Radome 

bulkhead 



070066 

Figure 1. Flight QS description. 


The Mock QS was designed with the purpose of replacing the flight QS in the build-up 
GVT approach. For the Mock QS to fulfill this purpose, it was necessary that its design meet the 
following requirements: 1) similar weight to flight QS weight, 2) similar CG to flight QS CG, and 
3) similar I xx for the third segment of the flight QS. A final design requirement imposed on the 
Mock QS was an interface with the Strongback and F-15B radar bulkhead in the same manner and 
with the same hardware as the flight QS installation assembly. A study was done that considered 
different materials, beam types, thicknesses, diameters, and lengths to determine the best design 
for the Mock QS. The Mock QS measures 19 ft and is shown in figure 2. 



13 ft 


6ft 


070067 

Figure 2. Mock QS description. 


7 


Several model updates coupled with a series of four GVTs were performed to acquire a 
correlated mated F 15B/QS FE model that could be used in flutter analyses. Figure 3 shows the 
various models of the Mock QS that were developed to accomplish this task. The first FE model 
created was a detailed FE model of the Mock QS, which was created in Patran and MSC.Nastran 
(MSC Software Corporation, Santa Ana, California) using the Mock QS computer aided design 
(CAD) models. It consisted of beam, plate, and solid elements and was fixed in a manner that 
represented how the flight QS mounted to the airplane. A modal analysis was performed on the 
detailed FE model. The results from this modal analysis rendered the target frequencies and mode 
shapes for the beam model, shown on the far right in figure 3. 


Mock QS CAD 


Mock QS detailed model 


Mock QS beam model 



Figure 




070068 


3. Mock QS model development. 


The original beam model was developed in Patran and MSC.Nastran. The lumped mass and 
beam elements were created manually based on the numerical mass properties of the detailed 
model. Because the model update process requires much iteration, using a reduced-order model 
decreases the computational effort required to match model structural properties to test data. 
This beam model decreased the computational effort in all analyses including the flutter analyses 
required at the conclusion of the model update process. Furthermore, matching the beam model to 
the detailed model first provided a good ‘initial guess’ for the Mock QS beam model and would, 
therefore, facilitate correlating the beam model to test data. 

The next step in this task was to correlate the Mock QS beam model to the actual Mock QS 
structure. This was accomplished in the first GYT, which consisted of attaching the Mock QS boom 
to a strongback and exciting it in the lateral and vertical directions as shown in figure 4. The results 
from this GVT were used to update only the structural properties of the beam model. The interface 
stiffness was not updated until the second GYT. Frequencies and mode shapes were extracted from 
all the following GVTs to provide target data for model validation (ref. 6). Attaching the Mock QS 
structure to the Strongback provided a nearly rigid foundation for the boom that ensured the boom 
modes were excited well. This also matched well with the boundary condition modeled. 
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Vertical shaker 


Mock QS boom 


Strongback 




Lateral shaker 




070069 

Figure 4. Mock QS ground vibration testing on the Strongback structure. 


When mating the updated Mock QS beam model to the updated airplane beam model, it was 
assumed that any discrepancies between the numerical and the test mode shapes and frequencies 
could be attributed to the unknown connection stiffness between the test article and airplane. The 
second GYT would empirically provide this connection stiffness. The airplane was constrained 
laterally and vertically, to achieve airplane modes that were decoupled from the Mock QS modes. 
The Mock QS was again excited both laterally and vertically. This GVT configuration is illustrated 
in figure 5. 
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Photo courtesy Tony Landis 

Figure 5. Mock QS boom ground vibration testing on the F-15B airplane. 


The third GYT consisted of the flight QS mounted on a strongback that provided the modal 
parameters to update the flight QS model. The fourth and final GVT served as a final check for the 
mated F-15B/QS frequencies and mode shapes. 

RESULTS AND DISCUSSION 

Because of the simplicity of the structure and the fidelity of the detailed FE model of the 
Mock QS, there was confidence in its representation of the actual Mock QS hardware. Therefore, 
in performing the first model update, the target modal values were provided by the numerical 
results of the detailed model to update the beam model. The first three modes of the detailed model 
or, target modes, are shown in figure 6. 


Mode 1 -lateral bending (5.675 Hz) Mode 2-vertical bending (6.602 Hz) 



Mode 3-torsion (7.876 Hz) 


070070 



Figure 6. Detailed Mock QS mode shapes and frequencies. 
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The first three modes of the detailed Mock QS were lateral bending at 5.68 Hz, vertical 
bending at 6.6 Hz, and torsion at 7.88 Hz. Table 2 quantifies the similarity in the modal parameters 
between the beam and detailed Mock QS models before the beam model was updated. 


Table 2. Comparison of frequencies and MAC values before update. 


Mode 

Detailed FE 
Model (Hz) 

Beam Model (Hz) 

MAC (%) 

Freq. Error (%) 

1 

5.675 

7.285 

(3 rd detailed mode) 

0.05 

28.37 

2 

6.602 

7.515 

(1 st detailed mode) 

0.02 

13.83 

3 

7.876 

8.977 

(2 nd detailed mode) 

0.03 

13.98 


Table 2 indicates that, initially, the beam model did not correlate well with the detailed model. 
The beam model modes were ordered differently: 1) torsion at 7.285 Hz, 2) lateral bending at 
7.515 Hz, and 3) vertical bending at 8.977 Hz. Consequently, this resulted in very low MAC values 
and high errors in frequency. The higher-than-expected frequencies imply that either the masses 
were underestimated or the stiffnesses were overestimated in the beam model. 

As discussed in the previous section, the first phase of updating the beam model consisted of 
updating the mass properties. The optimization variables were assigned upper and lower limits of 
±10 percent of the original values and included the mass values, CGs and inertias of the lumped 
mass elements. All mass properties were matched to within 0.5 percent error with the exception 
of the y-CG because the detailed model possesses a y-CG at zero. A comparison for all mass 
properties is shown in table 3. 


Table 3. Comparison of mass properties after update. 



Mass, 

lb 

x-CG, 

in 

y-CG, 

in 

z-CG, 

in 

hex’ 

lb-in 2 

h'Y’ 

lb-in 2 

^zz’ 

lb-in 2 

Kz’ 

lb-in 2 

^ZX’ 

lb-in 2 

Detailed 

400.01 

114.67 

0.00 

108.08 

2.60E+04 

2.16E+06 

2.15E+06 

7.40 

2.41E+04 

Beam 

400.59 

114.55 

0.03 

108.10 

2.60E+04 

2.16E+06 

2.15E+06 

7.38 

2.41E+04 

Error 

-0.14% 

0.10% 

NaN 

-0.02% 

-0.14% 

0.01% 

0.01% 

0.21% 

0.10% 
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The second phase of this model update consisted of producing an orthogonalized mass 
matrix. The variables were the same as in the previous phase; however, with the implementation 
of the constraint functions, the variables were limited to changing in such a way that the objective 
functions errors from the previous phase would not increase to a value greater than 8 (See table 1). 
At the conclusion of this phase, the off-diagonal term in the mass matrix of the greatest value was 
0.0135. The resulting orthogonalized mass matrix (in percent) was: 

100.00 -0.02 1.35 

-0.02 100.00 0.38 

1.35 0.38 100.00 


The third phase of implementing the DFRC mode matching code matches the beam model 
frequencies and mode shapes to the detailed model results. For this phase, the design variables 
were expanded to include element stiffness properties such as beam cross sectional area, area 
moment of inertias, and torsional constants. This set of design variables was also allowed to range 
±10 percent to complete the model update procedure. Table 4 demonstrates how the DFRC in- 
house mode matching code correlated the beam model to the detailed model. 


Table 4. Comparisons of frequencies after the first update. 


Mode 

Detailed FE 
Model (Hz) 

Beam Model 
(Hz) 

MAC (%) 

Freq. Error (%) 

1 

5.675 

5.674 

98.29 

-0.02 

2 

6.602 

6.607 

98.14 

0.08 

3 

7.876 

7.867 

78.64 

-0.11 


The updated beam model showed good correlation to the detailed model. The first two mode 
shape vectors were matched to less than two percent error, and the frequencies for all three modes 
were matched to within one percent error. The third mode shape, which was the torsion mode, 
did not possess as good a MAC value as the first two. Because of the geometry of the QS and the 
additional torsional stiffness provided by the radome, the first torsion mode of the actual QS does 
not appear until approximately 100 Hz. Therefore, it was deemed not critical to match the torsion 
mode with as much rigor as the lateral and bending modes, and a MAC value of 78 percent sufficed 
for this mode. 

Figure 7 shows the first three mode shapes of the Mock QS beam model after it was 
updated. 
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Mode 1 - Lateral bending, 5,674 Hz 



Mode 2 - Vertical bending, 6,607 Hz Mode 3 - Torsion, 7,667 Hz 



Figure 7. Mock QS beam model mode shapes and frequencies. 


The beam model was next updated to match test data that resulted from the first GYT. This 
first Mock QS GVT required that the Mock QS be mounted to a strongback. The first three modes 
of the Mock QS were 1) lateral bending at 5.583 Hz, 2) vertical bending at 7.752 Hz, and 3) torsion 
at 8.105 Hz. Updating the beam model to the detailed model first made updating the beam model 
to test data faster and less computationally intensive. Table 5 compares the beam model to the 
Strongback GYT results after only updating the mass properties. In this case, the only mass 
properties that were measured were total weight and x-CG. The numerical values were kept for all 
other mass properties. 


Table 5. Comparison of frequency and MAC values at update start. 


Mode 

GVT (Hz) 

Beam Model 
(Hz) 

MAC (%) 

Freq. Error (%) 

1 

5.583 

5.643 

99.29 

1.08 

2 

7.752 

7.529 

99.49 

-2.88 

3 

8.105 

8.142 

67.66 

0.45 


Because the beam model had already been through one update, its modal parameters already 
correlated well with the measured modal parameters. The design variables were treated similarly 
as before. They were assigned upper and lower limits of ± 10 percent of the original values and 
included the mass values, CGs and inertias of the lumped mass elements. The total weight and 
x-CG were updated to within a one percent error. 

The second phase of the model update, mass matrix orthogonalization, was also treated as 
previously done. The design variables were the same as in the previous phase. At the conclusion 
of this phase, the off-diagonal term in the orthogonalized mass matrix of the greatest value was 
0.0737. The resulting orthogonalized matrix (in percent) was: 
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100.00 

-1.54 

7.37 


-1.54 

100.00 

-6.17 


7.37 

-6.17 

100.00 


Finally, the beam model was updated to have frequencies and mode shapes that matched the 
measured data. The design variables were again expanded to include element stiffness properties 
such as beam cross sectional areas, area moment of inertias, and torsional constants that were 
allowed to range +10 percent. Table 6 shows how the DFRC in-house mode matching code 
correlated the beam model to the actual Mock QS boom. 


Table 6. Comparisons of frequencies after update to GVT data. 


Mode 

GVT (Hz) 

FE Model (Hz) 

MAC (%) 

Error (%) 

1 

5.583 

5.594 

99.25 

0.20 

2 

7.752 

7.764 

99.42 

0.15 

3 

8.105 

8.093 

88.51 

-0.15 


The resultant numerical modal parameters showed an excellent correlation to the measured 
modal parameters for the Mock QS test article. The first two mode shape vectors were matched 
to less than one percent error, and the frequencies for all three modes were matched to within 
0.5 percent error. Regarding the third mode shape, we see the same results as before. This mode, 
the torsion mode, did not possess as good a MAC value as the first two and was again, not regarded 
as being of great importance. 

The final model update required that only the spring element stiffnesses that connected the 
Mock QS hardware to the aircraft be updated. Initially, they were assumed to be nearly rigid 
and therefore, were assigned high values. The second GVT would provide their true values. The 
second GVT consisted of mating the Mock QS boom to the F-15B airplane. The airplane was 
constrained vertically by supporting it with aircraft jacks and laterally through use of an erector 
set. This configuration decoupled the aircraft modes from the boom modes. The measured modes 
and frequencies of this configuration were lateral at 6.06 Hz and vertical at 7.75 Hz. The torsion 
mode was not extracted for this configuration. Table 7 computes the frequency errors between the 
numerical and test modes. 


Table 7. Comparison of test and numerical modes of the mated Mock QS/F-15B airplane. 


Mode 

GVT (Hz) 

FE Model (Hz) 

Freq. Error (%) 

1 

6.064 

5.583 

8.62 

2 

7.756 

7.752 

0.05 
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The first airplane test mode was at a higher frequency than the first numerical mode by almost 
0.5 Hz. This was not expected, because the mated model already assumed a nearly rigid connection 
between the Mock QS and the airplane. The higher frequency of the first airplane test mode 
indicates that the airplane was stiffer laterally than the Strongback was, which can be explained 
by the high inertia of the airplane compared to the Strongback with the erector set providing a true 
lateral constraint. This also indicates that the assumption of a nearly rigid connection was accurate 
and is further supported by the fact that the second numerical mode frequency already matched the 
second test mode frequency to within a tenth of a percent. 

After the flight QS model was updated by test data, it was mated to the F-15B model using 
this nearly rigid connection stiffness. The project then proceeded to perform flutter analyses using 
the mated model. This final mated model matched well to GVT data, and therefore, no model 
update was necessary, and the flutter analysis did not need to be redone. The QS team could 
quickly proceed to the flight-testing portion of this project. 

CONCLUSIONS 

The Dryden Flight Research Center has developed an in-house mode matching code that 
updates both the mass and stiffness properties of a model. This method renders a more accurate 
model than the method of only updating stiffness properties. The three phases of optimization 
that the code employs are 1) updating mass properties, 2) orthogonalization of mass matrix, and 
3) updating frequency and mode shape. The code has been applied in the GAC QS project flown 
on the DFRC F-15B supersonic flight testbed. The build-up approach used to update the QS/F-15B 
mated model was completed using this code. During this project, it was used to perform two model 
updates. It successfully updated each model used for each mode of interest to within one percent 
error. This code was applied in such a way that allowed the project to proceed with flutter analyses, 
such that when the flight test article arrived at Dryden, the project could quickly move into flight 
testing. 
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